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ABSTRACT 



Aims. We present /, H, K spectrally dispersed interferometry with a spectral resolution of 35 for the Mira variable S Orionis. We aim 
at measuring the diameter variation as a function of wavelength that is expected due to molecular layers lying above the continuum- 
forming photosphere. Our final goal is a better understanding of the pulsating atmosphere and its role in the mass-loss process. 
Methods. Visibility data of S Ori were obtained at phase 0.78 with the VLTI/ AMBER instrument using the fringe tracker FINITO at 
29 spectral channels between 1 .29 fjm and 2.32 yum. Apparent uniform disk (UD) diameters were computed for each spectral channel. 
In addition, the visibility data were directly compared to predictions by recent self-excited dynamic model atmospheres. 
Results. S Ori shows significant variations in the visibility values as a function of spectral channel that can only be described by 
a clear variation in the apparent angular size with wavelength. The closure phase values are close to zero at all spectral channels, 
indicating the absence of asymmetric intensity features. The apparent UD angular diameter is smallest at about 1.3 /urn and 1.7 jim and 
increases by a factor of ~ 1.4 around 2.0yum. The minimum UD angular diameter at near-continuum wavelengths is UD =8.1±O.5 mas, 
corresponding to R ~ 420 R G . The S Ori visibility data and the apparent UD variations can be explained reasonably well by a dynamic 
atmosphere model that includes molecular layers, particularly water vapor and CO. The best-fitting photospheric angular diameter of 
the model atmosphere is Phot =8.3±O.2mas, consistent with the UD diameter measured at near-continuum wavelengths. 
Conclusions. The measured visibility and UD diameter variations with wavelength resemble and generally confirm the predictions by 
recent dynamic model atmospheres. These size variations with wavelength can be understood as the effects from water vapor and CO 
layers lying above the continuum-forming photosphere. The major remaining differences between observations and model prediction 
are very likely due to an imperfect match of the phase and cycle combination between observation and available models. 
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1 . Introduction ing of pulsation and mass loss. Observed radii of Mira stars have 

been found to differ for different optical and infrared bandpasses 
Mira stars are low-mass, large-amplitude, long-period variable ( _ Thompson et aL [JooH Mennesson et al. 2002; Ireland et 
stars on the AGB, evolving toward the planetary nebula and a , m ^ Perrin e{ a , g^. Eifjner e{ aJ 2QQ7h and Mf , hafj 



d ; white dwarf phases. They exhibit a mass-loss rate on the order beefl attributed to the ce of moleC ular layers located above 

of~10- 6 M G /yearthatsigmficantlyaffectsthefurtherstellarevo- ±e continuum .f orming photosphere. Here, we present both a 

lution and is one of the most important sources for the chemical spectro -interferometric observation of the Mira star S Ori that 

enrichment of the interstellar medium. The dust condensation se- coyers ^ near . infrared 7) ^ and K bands simu i tane ously at a 

quence, the wind-driving mechanism, and the role of pulsation ^ resolution of 35 and a compal - ison t0 rece nt self-excited 

are currently not well understood in particular for o = ch d ic d j at heres . 

AGB stars (Woitke et al. I2006I Hofner & Andersen 2007). The 3 F 
pulsating atmospheres of Mira stars can become very extended 

because of dynamic effects including shock fronts, and they are S Ori is a Mira variable star with spectral type M6.5e-M9.5e 

very cool in their outer parts. Here, molecules can form, which and v magnitude 7.2-14.0 (Samus et al.|2004). We use a Julian 

for O-rich stars are most importantly H 2 0, CO, TiO, and SiO Date of last maximum brightness T = 2453190 days, a period 

(Ts uji eta l. [19971 Tej et al. 120051 Ohnaka 12004b . Wittkowski et p = 430 da y s ( as in Wittkowski et al. |2007|l, and the distance 

al. (120071) found for the case of S Ori that A1 2 3 dust condenses of 480 P c ± 120 P c from van Belle et al - <2002). The broad- 

within the extended atmosphere at phase-dependent distances of band near-infrared K UD angular diameter of S Ori h as been 

1.8-2.4 photospheric radii. This extended atmosphere, which is measured by van Belle et al. ( 1996), Millan-Gabet et al. (|2005 ), 

characterized by phase-dependent temperature and density strat- and Boboltz & Wittkowski (12005|l to values between 9.6 mas and 

ifications, the presence of molecular layers, and the formation 10.5 mas at different phases. Joint VLTI /MIDI and VLBA/SiO 

of dust, is thus of particular interest for our better understand- maser observations by Wittkowski et al. Q2007|l have shown that 

S Ori exhibits significant phase-dependencies of the atmospheric 

* Based on observations made with the VLT Interferometer (VLTI) extension and dust shell parameters with photospheric angular 

at Paranal Observatory under program ID 080.D-0691 diameters between 7.9 mas and 9.7 mas. 
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Table 1. Observation log. Night starting 12 October 2007, JD 2454386. 



Target 


Purpose 


©LD 


DIT 


Time 


®Vi s 


B„ [m] 


PA„ 


AM 


Seeing 


To 






[mas] 


[msec] 


[UTC] 




E0-G0/G0-H0/E0-H0 


deg 




["] 


[msec] 


45 Eri 


Calibrator (K3 II-III) 


2.15±0.04 


25 


08:03-08:07 




16.0/31.9/47.9 


-107 


1.1 


1.3 


1.3 


45 Eri 


Calibrator (K3 II-III) 


2.15±0.04 


50 


08:09-08:12 




16.0/32.0/48.0 


-107 


1.1 


1.3 


1.3 


y Eri 


Check star (M0.5 Illb) 


8.74±0.09 


25 


08:22-08:26 




15.6/31.2/46.8 


-104 


1.1 


1.4 


1.2 


7 Eri 


Check star (M0.5 Illb) 


8.74±0.09 


50 


08:28-08:32 




15.5/31.0/46.5 


-104 


1.1 


1.4 


1.2 


SOri 


Science target 




25 


08:45-08:49 


0.78 


15.9/31.8/47.8 


-107 


1.1 


1.3 


1.3 


SOri 


Science target 




50 


08:52-08:57 


0.78 


16.0/31.9/47.9 


-107 


1.1 


1.3 


1.3 


a Hor 


Calibrator (K2 III) 


2.76+0.03 


25 


09:14-09:18 




14.5/28.9/43.3 


-91 


1.1 


1.2 


1.4 


a Hor 


Calibrator (K2 III) 


2.76±0.03 


50 


09:21-09-24 




14.3/28.7/43.0 


-90 


1.1 


1.4 


1.2 


SOri 


Science target 




25 


09:36-09:40 


0.78 


15.9/31.8/47.7 


-106 


1.1 


1.5 


1.1 



2. Observations and data reduction 

We obtained near-infrared J, H, K interferometry of S Ori with 
the instrument AMBER (Petrov et al. 120071 ) in low-resolution 
mode at the ESO VLTI using the fringe tracker FINITO and 
three VLTI Auxiliary Telescopes (ATs) on 12 October 2007 (JD 
2454386). The ATs were positioned on stations E0, GO, and HO. 
The date of observation corresponds to a visual phase cDvis=0.78, 
with an uncertainty of about 0.1. The details of the observing 
sequence are listed in Table [T] including the projected baseline 
lengths (B p ) and position angles (PA ;) ), the airmass (AM), and 
the optical seeing and coherence time. Ambient conditions were 
not very good but stable. The airmass was the same for all ob- 
servations. Data were recorded using two different detector inte- 
gration times (DITs) of 25 msec and 50 msec. 

The VLTI fringe-tracker FINITO records fringes on the two 
shortest baselines using 70% of the //-band light. Output signals 
are processed in realtime and used to compensate for the fringe 
motion due to atmospheric turbulence. Owing to the low coher- 
ence time during the observations, FINITO was only able to pro- 
vide average performances of 0.2-0.4yum RMS (to be compared 
with 0.1 /im RMS achieved in good conditions). Nevertheless 
such a performance is sufficient for increasing the signal-to- 
noise ratio (S/R) of the AMBER data and stabilizing the transfer 
function. All observations used here were obtained with exactly 
the same FINITO controller parameters, which is important for 
avoiding systematic biases in calibrating the absolute visibility. 

In addition to S Ori, the two calibration stars 45 Eri and 
a Hor were observed close in time, as was y Eri. The last is a 
regular non-pulsating M 0.5 giant with a well-known angular di- 
ameter similar to that of S Ori which is not expected to exhibit 
strong effects from molecular layers. This data is used to check 
that any strong wavelength-dependent features found for S Ori 
are not caused by any systematic effects of the instrument. The 
spectral types and angular diameters of the calibration and check 
stars are from Borde et al. (2002). 

Raw visibility and closure phase values were computed us- 
ing the latest version of the amdlib package (version 2.0 beta 
2b) and the yorick interface provided by the AMBER consor- 
tium and the Jean-Marie Mariotti Center. Absolute wavelength 
correlation was performed by correlating the raw spectra of all 
four stars with a model of the atmospheric transmission with the 
same spectral resolution. In particular, we used a plateau in the 
transmission curve at A ~ 2.0 /mi. The offset with respect to the 
original wavelength table was 3 spectral channels (0. 1 jum at a 
wavelength of 2.0/mi). We estimated the error of the absolute 
wavelength calibration to 1 pixel (~ 0.03 yum). Individual frames 
were averaged after frame selection keeping 70% of the best 
frames based on piston (to remove the frames when the FINITO 
loop was not closed) and out of these keeping 30% of the best 



frames based on S/R. We verified that keeping up to 80% of the 
best frames based on SNR did not significantly change the re- 
sults. The S Ori and y Eri visibility data were calibrated for each 
DIT value separately using the 45 Eri and a Hor calibration star 
data. After calibration, the different calibrated S Ori and y Eri 
data were averaged. The errors of the calibrated visibility data 
include the statistical error of averaging the single frames, the 
errors of the calibration stars' angular diameters, and the varia- 
tion of the available transfer function measurements. 



3. Results 

Figures [TJ and [2] show the resulting visibility and closure phase 
data for S Ori and for the check star y Eri, respectively. The gaps 
in the visibility data around 1.45 //m and 1.85 fim correspond to 
the regions between the bands. Also shown are the best fitting 
models of a UD with a constant diameter, of a Gaussian disk 
with a constant diameter, and of atmosphere models. The lat- 
ter are for S Ori the M18n model (described in detail below in 
Sect. |4), and for y Eri an ATLAS 9 model atmosphere (Kurucz 
119931 ) with r e ff=3750K, log g=1.5, and solar chemical abun- 
dance as expected for y Eri (Borde et al. 120021 ). The comparison 
of the S Ori visibility data to the dynamic model atmosphere is 
described in detail below in Sect. [4] 

The calculation of synthetic visibility values and the fits to 
the interferometric data were performed as in Wittkowski et al. 
(2006, 2007 ). The fitted angular diameter of the Gaussian model 
corresponds to the FWHM, that of the plane-parallel ATLAS 9 
model to the 0% intensity (limb-darkened) radius, and that of the 
M18n model to the 1.04/vm (photospheric) radius (as defined in 
Ireland et al. 2004b). The best-fit angular diameters are 
S Ori: UD = 10.8 mas; G au.ssian=6.9mas; M i8n=8.3mas 
y Eri: ©uD=8.5mas; ©ATLAS9=8.9mas. 
The errors are cr ~ 0.2 mas 

The visibility data of the check star y Eri can be described 
well by a UD of constant diameter and by the ATLAS 9 model 
atmosphere. There are no significant wavelength-dependent de- 
viations between measured visibility data and the UD model. It 
is not yet clear whether the relatively low /-band visibilities for 
the G0-H0 baseline and the deviations in the closure phase val- 
ues near the flip are caused by an asymmetric stellar surface fea- 
ture or a systematic calibration uncertainty. Note that the angular 
diameter of y Eri based on a given model is well-constrained by 
the position of the visibility minimum, which is independent of 
an absolute visibility calibration. The resulting angular diame- 
ter ©ATLAS9=8.9±0.2mas is consistent with the value given in 
Borde et al. of ©ld=8.74±0.09 mas. 

The visibility data of S Ori show significant wavelength- 
dependent features clearly deviating from UD and Gaussian 
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Fig. 1. Measured S Ori visibility data compared to models of a UD with a constant diameter (red dashed lines), of a Gaussian disk 
of constant diameter (green dashed line), and of the M18n atmosphere model (blue solid line). For the projected baseline lengths 
and angles see Table Q] 
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Fig. 2. Measured y Eri visibility data compared to models of a UD with a constant diameter (red dashed line) and of an ATLAS 9 
model atmosphere with r e ff=3750K, log §=1.5, solar chemical abundance (blue solid line). For the projected baseline lengths and 
angles see Table Q] 
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Fig. 3. S Ori UD diameter values as a function of wavelength 
compared to the prediction by the M18n model atmosphere. Also 
indicated are the positions of H2O and CO bands after Lancon 
& Wood (120001 



models of constant diameter on all three baselines. This indi- 
cates variations in the apparent angular diameter. The S Ori clo- 
sure phase values are consistent with zero at all spectral chan- 
nels, indicating the absence of asymmetric features in the inten- 
sity distribution. In order to characterize the variation of S Ori's 
angular diameter as a function of wavelength, we fitted UD di- 
ameters to the data of each spectral channel separately. Figure [3] 
shows the resulting UD diameter values as a function of wave- 
length. Note that the intensity profile of a Mira star is generally 
not expected to be close to a UD and that this approach can only 
give a rough estimate of S Ori's characteristic size as a function 
of wavelength. Fits of Gaussian functions lead to similar results. 
The apparent UD angular diameter shows clear variations with 
wavelength. It is smallest at about 1.3//m and 1.7 /mi and in- 
creases by a factor of ~ 1.4 around 2.0//m. The minimum UD 
angular diameter of S Ori at the near-continuum wavelengths is 
©ud=8. 1+0.5 mas. This result at visual phase 0.78 is consistent 



with the S Ori photospheric angular diameters between 7.9 mas 
at phase 0.55 and 9.7 mas at phase 1.15 derived in Wittkowski 
et al. (120071 1 based on VLTI/MIDI data and modeling with dy- 
namic model atmospheres and a radiative transfer model of the 
dust shell. These photospheric angular diameters are also con- 
sistent with previous broadband measurements (cf. Wittkowski 
et al. 120071) . With the adopted distance to S Ori, the angular di- 
ameter ©ud=8. 1+0.5 mas corresponds to aradius418+ 130 R Q . 

4. Comparison to dynamic Mira star atmosphere 
models 

Few dynamic atmosphere models for oxygen-rich Mira stars are 
available that include the effects of molecular layers. The P and 
M model series (Ireland et al. 2004a 2004b) are complete self- 
excited dynamic model atmospheres of Mira stars designed to 
match the prototype oxygen-rich Mira stars o Cet and R Leo. 
They have been used successfully for comparisons to recent 
broadband interferometric data of o Cet and R Leo (Woodruff et 
al. 120041 Fedele et al. 120051 Compared to o Cet and R Leo, S Ori 
is a slightly cooler Mira variable with a longer period, a higher 
main-sequence precursor mass, and a larger radius. However, 
when scaled to variability phases between and 1 and to the 
corresponding angular size on the sky, the general model results 
are not expected to be dramatically different for S Ori compared 
to o Cet and R Leo (cf. the discussion in Wittkowski et al. 120071 1. 
The M model series was chosen to model the atmosphere of 
S Ori by Wittkowski et al. (2007) as the currently best avail- 
able option to describe Mira star atmospheres. Monochromatic 
center-to-limb variations (CLVs) at 46 angles between and 5 
R p based on the P and M models were recomputed for the wave- 
length range from 1-2.5 fim in steps of 0.001 /mi. 

The P and M dynamic model atmospheres predict signifi- 
cant changes in the monochromatic radius Rj = Ra(ta - 1), 
caused by molecular layers that lie above the continuum-forming 
photosphere and significantly affect certain bandpasses. Figure|4] 
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Fig. 4. Monochromatic T\ = 1 Mira star radii predicted by the M model series for the example of the M16n (model phase 1.60), 
M18n (1.84), and M20 (2.05) models. The red line indicates the model continuum radius excluding atomic and molecular features. 
The spectral resolution of the model is 0.001 //m. Also indicated are the positions of H2O and CO bands after Lancon & Wood 
(2000). 



shows the monochromatic radius Rx = Ra(?a = 1) in units of the 
non-pulsating parent star radius for the example of three mod- 
els of the series. It illustrates the strong phase dependence of 
the molecular layers. The models are M16n (model phase 1.60), 
M18n (1.84), and M20 (2.05) models. The red lines denote for 
comparison the Ra values solely based on the continuum radi- 
ation excluding all atomic and molecular features. Also shown 
are the positions of the H2O and CO bands after Langon & Wood 
(2000) and references therein. The most prominent features of 
these model curves in the near-infrared region are two water va- 
por features around 1.4//m and 1.9//m, and also CO features 
around 1.6//m and 2.4 /mi. The strengths, shapes, and widths 
of these molecular features depend strongly on the stellar phase 
(and also on cycle), as is evident from the comparison of the 
three model curves. Also, the relative strengths of the molecular 
features varies with stellar phase. 

Model M18n provides the best formal fit to our S Ori 
AMBER visibility data out of the 20 available phase and cy- 
cle combinations of the M series. The synthetic visibility values 
based on the M18n model compared to our AMBER observation 
are indicated in Fig.Q] Here, the angular photospheric diameter 
corresponding to the 1 .04 /mi (photospheric) model radius (de- 
fined in Ireland et al. 12004b) is ©Mi8n=8.3 + 0.2 mas, consistent 
with the UD diameter at near-continuum wavelengths 1.3 /mi 
and 1.7 fim of 0ud=8. 1+0.5 mas. The theoretical Ra(ta = 1) 
radii in Fig. |4]cannot be compared directly to the UD diameters 
derived from our AMBER data (Fig. [3]), because of the different 
spectral resolution and because the model-predicted CLVs can 
be very different from a UD model (so that different radius def- 
initions are not equal). The translation of the model prediction 
into a UD diameter depends on the exact shape of the bandpass- 
averaged CLV and the baselines used. To compare the model 
predictions to the measured UD values in Fig. [3] we fitted UD 
diameters to the synthetic visibility values of the M18n model 
using exactly the same spectral channels, baseline configuration, 
and fit method as for our AMBER data. The resulting model 
prediction for the UD diameter as a function of wavelength is 
shown in Fig. [3] in comparison to the values derived from the 
observation. 

Figures Q] and [3] show that our AMBER visibility data can 
be described reasonably well by the dynamic atmosphere model 
M18n. The differences between observations and model pre- 
diction are most likely due to an imperfect match of the phase 
and cycle combination between observation and available M 
models of the series. Looking at Figs. [3] and [4] it is evident that 
a better fit to our AMBER data could be obtained with a model 
that shows a stronger 1.4 /tm water vapor feature of the same 



shape compared to M18n (as seen for instance in the case of 
M20), and at the same time a just as strong but broader 1.9/mi 
feature (as for instance in the case of M16n). It is quite possible 
that such a combination of the two water-vapor features could 
appear for a model of another phase-cycle combination. Also, 
some differences between M model predictions and observations 
of S Ori are expected due to the different stellar parameters of 
S Ori compared to the parent star of the M model series. Finally, 
differences can also be caused by remaining uncertainties in 
the absolute calibration of visibility values and of the wave- 
length scale. The S Ori dust shell of ry=1.5-2.5 as modeled in 
Wittkowski et al. (2007) is not expected to have a noticeable 
effect on our near-infrared visibility data, because its contribu- 
tion to the visibility was already small at 8 //m and because the 
visibility data are calibrated separately for each spectral channel. 

In summary, our AMBER observations of S Ori generally 
confirm the predictions by the M model series and we find that 
the observed variation of diameter with wavelength can be un- 
derstood as the effect of phase-dependent water vapor and CO 
layers lying above the photosphere. The M model series can 
be used reasonably well to model the atmosphere of a Mira 
star such as S Ori and to derive a reliable photospheric radius 
based on broadband data. More such observations are needed to 
confirm and constrain the model predictions in more detail and 
to monitor the predicted phase dependence of the strength and 
characteristics of the molecular layers. Simultaneously obtained 
spectra would be a valuable addition. 

Acknowledgements. We acknowledge with thanks the use of the AMBER data 
reduction software from the JMMC (version 2.0 beta 2b). 
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